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ABSTRACT 

Recent simulations have indicated that the dark matter halos of galaxy clusters should feature 
steep density jumps near the virial radius. Since the member galaxies are expected to follow similar 
collisionless dynamics as the dark matter, the galaxy density profile should show such a feature as 
well. We examine the potential of current datasets to test this prediction by selecting cluster members 
for a sample of 56 low-redshift (O.I < z < 0.3) galaxy clusters, constructing their projected number 
density profiles, and fitting them with two profiles, one with a steep density jump and one without. 
Additionally, we investigate the presence of a jump using a non-parametric spline approach. We find 
that some of these clusters show strong evidence for a model with a density jump. We discuss avenues 
for further analysis of the density jump with future datasets. 


1. INTRODUCTION 

Galaxy clusters contain a representative sample of the 
matter in the universe: the dominant constituent is dark 
matter, while the baryonic components incl ude hot ga s 
and the galaxies themselves (for a review, see lVoitir2nn5h . 
A longstanding analytical prediction for cosmological 
structure formation is the existence of a shock bound¬ 
ing the hot, gaseous intracluster medium and a corre¬ 
sponding jump in the dark matter profile, co incident with 
the v irial radius of the galaxy cluste r fe.g.. iBertschingei 
[198^. More recen tly, sim ulations bv iDiemer fc Kravtsov 
(|2014ll (hereafter iDKlJf l have reinforced the expecta¬ 
tion that the dark matter halo itself should exhibit a 
sharp density steepening near the viri al radius (see also 
[Adhikari. Da.lal. fc Chainberlainll2014fl . 

The member galaxies of a cluster are expected to trace 
the cluster’s dark matter profile in the cluster outskirts, 
as they are subject to similar collisionless dynamics; ac¬ 
cordingly, we may expect to see such a feature not just 
in the dark matte r prof ile but also in the galaxy den¬ 
sity profile. iTullvl (l2010f) examined the distributions of 
the galaxies in the Coma and Virgo clusters and de¬ 
tected sharp density cut-offs at radii of 3 and 2 Mpc, 
respectively, from the cluster centers, which he identified 
with the caustics of second turnaround. These corre- 
spond roughly to the virial radii of these c lusters (e.g., 
iKubo et al.ll2007HKarachentsey et aI1l2014ll . The analy- 
sis of iTrentham &: Tullvl (|2009ll identified a similar fea- 
t ure in the galaxy group NGC 1023. 

More. Diemer. fc Kravts^ (|20I5[ 1 mentioned that 
hints of such a jump ma y hav e been seen in other 
data sets. For instance, (IMH) measured the 

second turnaround radius, which is associated with 
the aforementioned density jumps, of s e veral groups 
and cluste rs. Additionally, iRines etldl (j20I3l l (here¬ 
after lR,I3f l presented spectroscopic velocities for clus¬ 
ter members in 58 galaxy clusters, from which they re¬ 
constructed the density profiles, which appear to show 


a deficit of galaxies with resp ect to an NEW profile 
(jNavarro. Frenk. fc Whit3Il997fl . However, they noted 
that their results at large radii from the cluster center 
may be affected by having few spectroscopically observed 
cluster members in the exterior re gions. 

In this paper, we focus on the IRI3I sample of clus¬ 
ters and use public data, from the Sloan Digital Sky Sur¬ 
vey (SDSS) ([York et al.ll2000[ l to select cluster members 
photometrically. Rather than use spectroscopic veloci¬ 
ties, we construct the radial number density profiles of 
member galaxies and fit them with two functional forms 
to examine whether we can detect a feature consistent 
with the predicted density jump in any of the clusters 
using available data sets. Where applicable, we assume 
the standard AGDM cosmology with Ha = 0.73 and 
p,m. = 0.27, consistent with the parameters selected by 
IDK14I to enable direct comparison. To compare results 
across clusters we also use the measure Ra , which is the 
radius within which the mean mass density p = Api,{z), 
where pb is a specified cosmological background density 
and A is the density contrast with respect to Pb{z). 

2. DATA 

2.1. SDSS Catalogs 

To prob e the existence of the density jump predicted 
bv lDKl4 , we focus our attention on the sam ple o f 56 low- 
redshift (0.1 < z < 0.3) galaxy clusters from lRlSil (of their 
original 58 clusters, we do not include MS0906/A750 in 
our samples, since these clusters are nearly coincident 
and hence it is difficult to make a clean selection of mem¬ 
ber galaxies for each). Since their cluster sample was se¬ 
lected from regions of sky covered by SDSS, we use pub- 
licl y available data fr om the most recent Data Release 
12 (|Alam. et al.ll201^ to select cluster members. 

We obtain photometric catalogs from the SDSS Sky- 
Server SQL servei0. For each cluster, we query galax¬ 
ies within 1.5 degrees of the cluster center (which, 
from IR13I . is the X-ray center) in both right ascension 
(RA) and declination. We further restrict our query 
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each galaxy, we obtain the dereddened ‘model’ mag¬ 
nitudes (which we will use throughout the remainder 
of thi s work’), corrected for Galactic extin ction accord¬ 
ing to iSchlegel. Finkbeiner. fc DavisI (jl998ll by the SDSS 
pipeline, and we also select photometric redshift esti¬ 
mates and associated redshift quality estimates from the 
‘Photo-z’ table. 

2.2. Cluster Member Selection 

We will test two cluster member selections. The first, 
which we will refer to as Selection A, is based solely on 
photometry, primarily comprising red cluster galaxies se¬ 
lected from their location on a color-magnitude diagram 
but also including some galaxies with sufficiently secure 
photometric redshift estimates that are blueward of the 
red sequence. The second, which we call Selection B, 
folds in the spectroscopy of iRl.'ll to include additional 
cluster members and to reject some of the photomet¬ 
rically selected galaxies whose redshift estimates place 
them beyond the cluster. These selection procedures are 
outlined in more detail in the following sections. 

2 . 2 . 1 . Red Sequence Cluster Member Selection 

Both of our methods of cluster member select ion are 
based upon the target selection proces s of IrTI which 
relied on the red sequence method of iGladders fc Yed 
(j2000f) . The red sequence refers to a linear feature in 
the color-magnitude diagram of galaxies in a cluster field 
that arises from the population of red, early-type galax¬ 
ies that have been observed to comprise the majority of 
cluster members. By selecting galaxies within a range 
around this l ine, w e can obtain a sample of cluster mem¬ 
bers. As in mi we construct a diagram of g — r vs 
r and fix the slope of the red sequence line to —0.04. 
We then select the appropriate intercept by examining 
the distribution of galaxies within 5 arcminutes of the 
cluster center. We can also use the spectroscopic data 
publicly available from IR13I : by matching their spectro¬ 
scopic catalogs to SDSS photometry, we can additionally 
calibrate our red sequence selection, as will be described 
in Section [2.2.31 

Having thus identified the linear red sequence feature, 
we select galaxies that are within an offset of 0.1 magni¬ 
tude above the red sequence line, and 0.15 magnitude be¬ 
low the red sequence line, as illustrated in Figure[T] This 
is in contrast to the choice of ±0.3 magnitude in IR13L 
which was found to be much wider than the actual red 
sequence as indicated by the spectroscopy. The asym¬ 
metry of the limits above and below the red sequence is 
moti vated by the results of Section 5.2 and Figure 16 of 
IR13I . which suggest that there are fewer cluster members 
above the red sequence than below it, as confirmed by 
the righthand panel of Figure [l] which shows the distri¬ 
bution of spectroscopically confirmed members. 

2 . 2 . 2 . Contribution of Non-Red Sequence Galaxies 

Since the red sequence method selects a specific pop¬ 
ulation of galaxies, it is worth examining whether the 
exclusion of other types of galaxies i n the cluster affects 
our result. The analysis of lDressl^ (|1980f) . which exam¬ 
ined the galactic populations of over 50 galaxy clusters, 
found that blue, late-type galaxies comprise a small frac¬ 
tion of the populations of rich galaxy clusters, although 


this fraction increases as the density of galaxies decreases. 
This is supported by the spectroscopic data provided by 
El, who find that within the virial radius of the clus¬ 
ter the fraction of blue galaxies is < 10%, while over the 
entire radial range for which they obtained spectra, the 
fraction is < 15%. 

Since we have photometric redshift estimates, we use 
these to also select cluster members that are blueward of 
the red sequence. We restrict our attention to galaxies 
that have secure photometric redshifts, as indicated by 
the ‘photoErrorClass’ and ‘zErr’ tags; the former we se¬ 
lect to be equal to 1, and the latter we restrict to be less 
than 0.03. We then mark as additional cluster members 
galaxies that were not originally selected via the red se¬ 
quence method and whose photometric redshift e stima te 
is within 0.03 of the cluster redshift provided bv lR,13l . 

The population of galaxies with photometric redshift 
estimates of sufficient quality as indicated by the SDSS 
pipeline is too small for most clusters to use as the basis 
for our analysis; within the central 1.5 /i“^Mpc, there 
are on average only about 60 galaxies with photometric 
redshifts satisfying \zc — Zg\ < 0.03 and subject to the 
above quality cuts, many of which were already picked 
up via the red sequence method. Additionally, basing a 
selection on these redshifts would exclude most cluster 
members with r > 19. Accordingly, as noted above, we 
must combine them with the galaxies selected via the red 
sequence method. This final selection, denoted Selection 
A, for an example cluster is shown in the lefthand panel 
of Figured 

2 . 2 . 3 . Additional Tests of the Cluster Member Selection 

While neithe r the SDSS photometric redshifts nor the 
spectroscopy of lR13l identify enough cluster members for 
our analysis, we can use this additional information to 
refine the red sequence selection in a second selection. 
As before, after selecting red sequence cluster members 
using the procedure outlined in Section 12.2.11 we add in 
bluer cluster members based o n SDS S photometric red¬ 
shifts as described in Section 12.2.21 However, in this 
case, we can further refine the selection by also reject¬ 
ing cluster members if the SDSS photometric redshifts 
are sufficiently secure and the members in question have 
redshifts in the range |zc ~ Zg\ < 0.03, where Zc is the 
cluster redshift and Zg is the photometric redshift esti¬ 
mate. 

We can use the spectroscopic data of lR,13l similarly. Af¬ 
ter matching their tables to the SDSS catalogs, we add 
in cluster members as identified by the spectroscopy and 
reject galaxies that were originally selected by the red 
sequence method or phot ometr ic redshifts but are iden¬ 
tified as non-members by IR,13L This method of cluster 
member selection, which we term ‘Selection B,’ is sum¬ 
marized in the righthand panel of Figure [2l This selec¬ 
tion is more observationally expensive than Selection A, 
so it is worth testing both methods to see whether the 
additional information makes a difference in the results. 
On average, the numbers of cluster members selected by 
these two methods differ by about 5% within 1.5 /i“^Mpc 
of the cluster center. We will test the effect of this selec¬ 
tion on our analysis in Sectional 

As an additional step in verifying our cluster mem¬ 
ber selection, we construct a projected radial density 
plot consisting of the total radial density profiles of all 
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Fig. 1.— The red sequence cluster member selection for example cluster A1068. Left: in the color-magnitude diagram, we initially plot 
the location of galaxies within 5 arcminutes of the cluster center, which should be dominated by cluster members, to pr ovide a n indication 
of the red sequence. Middle: The galaxies selected by the red sequence method using the limits described in Section 12.2.11 we plot the 
galaxies thus selected that are within 2.5 h~^M.pc of the cluster cent er. The labelled intercept (g — r)i is selected to be the value at r = 16. 
Right: The cluster members as selected via the spectroscopy of IR13I. which can be used to verify the red sequence selection. 


sources in our catalog (subject to r < 20 ) as well as 
the non-cluster density profile, computed by subtracting 
the counts of cluster members from the total number of 
galaxies in each bin. In the total profile, we expect to 
see a steep increase in galaxies at small radial distances 
that flattens out at large R. After subtracting out the 
contribution from these cluster members, the resulting 
non-cluster member profile should be roughly flat. Fig¬ 
ure [3] shows these density plots for four example clusters, 
using Selection A. 


3. METHOD 
3.1. Approach 


Our fiducial analysis relies on cluster members selected 
as described in Section [2 subject to a magnitude cut of 
r < 20. We determine the number counts of galaxies in 
bins of 0.1 /i“^Mpc, from which we construct the pro¬ 
jected number density profile of member galaxies, N{R), 
with R indicatii ig the projected radius, using the cluster 
centers listed in lR13l . We use Poissonian error bars. 

We employ two methods to test for evidence of the 
density jump. In the first, we fit t he profiles using the 
fitting formula provided by IDK14I to incorporate their 
predicted density steepening: 


noK^r) = nin(r) 



-y 


5i?200 


-f 1 


( 1 ) 


As suggested bv lDKl4 we fix /3 = 6 and 7 = 4, a choice 
that yields a dependence of on the mass accretion rate, 
F: 


n = (o.62 -h 1.186-2^/3^ X i?200 (2) 


We further select the NFW profile with two parameters, 
ris and Cs, as our inner density profile riin: 


?t-nfw(^) 


__ 

t/ts {l + r/rsf 


(3) 


We n ote that IDK14I used the Einasto function (lEinastd 
1196511 instead of the NFW for the inner profile, but in 
the regime of interest, the distinction between the two is 
negligible. We will refer to the model given by Equation[T] 
with nin(r) = nNFw(' 6 ) as the ‘Density Jump’ model, or 
‘DJ’ model in abbreviation. 

Since the density jump feature occurs around the virial 
radius, we use a fitting range of i? < 2 i? 200 : which is 
beyond the range that is typically fitted well by an NFW 
profile (we als o exclud e the inner R < 0.1i?vir in the fit, 
consistent with lDKT4l . Accordingly, the two models that 
we will compare via fitting will be the DJ profile and a 
profile given by an NFW profile with an outer term: 


n(r) = nyiFw(r) + 


be 


5i?200 


-f 1 


(4) 


These profiles are projected numerically and fitted to 
our galaxy density data. The free parameters in our 
fit are rig, r^, be, and Sg for both profiles; the full DJ 
formula has one additional free parameter, F. We re¬ 
strict the fits to the range 0 < F < 5 and to rea¬ 
sonable values of rig and rg, the latter of which is re¬ 
stricted by the value of the NFW concentration pa¬ 
rameter, c = i? 200 c/Fs. Observations have suggested 
that this val ue is lower for galaxy profiles than for dark 
matt e r le.g..lLin. Mohr fc Stanford 12004 iHansen et al.l 
120051 iBudzvnski et al.l 1201211 . and we set the range as 
2.0 < c < 6.0 in our fits. We refer to the Appendix of 
IDK14I for appropriate ranges for Se (0.5 — 2.0) and be 
(0.1 — 4.0). We additionally fix Um, which will be dis¬ 
cussed further in Section [321 and we select as ou r uppe r 
limit of integration R — 10i?vir, same as used bv lDKl4 
We then examine the results of the fits using the 
Akaike Information Criterion (AIC; lAkaikel 11974^ an d 
the Bayesian Information Criterion fBIC: ISchwarzlll978ll , 
which provide a means of comparing models fitted to 
data. As the use of these methods in astrophysics and 
cosm ology has been discussed in a number of papers 
Ce.g.. lTakeuchill200(A iLiddldl^OOTt [Broderick et al.ll2011l : 
[Tan fc Biswasll2012f) . we simply mention the most salient 
qualities here and refer the reader to these works for fur- 
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Fig. 2. — Top: The cluster member selection for example cluster 
A1068 using Selection A, which comprises galaxies selected using 
the red sequence method and galaxies with SDSS photometric red- 
shifts, which add in some galaxies blueward of the red sequence. 
Bottom: The cluster member selection for this same cluster using 
Selection B, which combines the red sequence, photometric red- 
shifts, and spectroscopic redshifts. 

ther details. To apply these criteria, we compute the 
following statistics for each fit: 


AIC = x^+2p + 


+ 1 ) 

A/” — p — 1 ’ 


BlC = x^+ pH-^), 


(5) 

( 6 ) 


where p is the number of parameters in the fit, Af is 
the number of data points being fitted, and x^ is the 
standard minimized goodness-of-fit parameter. In the 
case of the AIC, we have also included a correction term 
of 2 p(p + 1 )/ ( Af — p — 1) which is recommended fo r small 
values of A/” (|Burnham fc AndersonI I2002L |2004|1 . The 


model that is preferred by these criteria is the one with 
the lower IC = AIC, BIC value. If we compute AIC = 
IChigh — ICiow, then, roughly, values of AIC = 1 — 5 are 
indicative of ‘positive’ evidence in favor of the model with 
lower IC and values of AIC > 5 denote ‘ strong’ evidence 
(e.g., ILiddldl^OOTt IBroderick et al.l[2011ll . 

Our second method consists of smoothing the profiles 
by fitting a smoothing cubic spline to our data. For a 
data set with measured values pi and errors ai at a set 
of points ri, the smoothing spline /(r^) is constructed to 
satisfy the condition 



where S' is a constant that interpolates between smooth¬ 
ing and fitting: that is, when S = 0, the spline is 
forced to pass through every data point, so that there 
is no smoothing, whereas as S is increase d, the curve be ¬ 
comes smoothe r at the expense of the fit (jde BooHl200lD . 
iReinschl (|1967ll argues that the smoothing parameter S 
should be chosen in the range A/”—•\/2A/’ < S < Af+V^Af, 
where Af is the number of data points over which we 
construct the spline, if the at are estimates of the stan¬ 
dard deviation in pi. We accordingly choose three val¬ 
ues within this range to compare to an NFW fit: S = 
Af — V2Af, Af, Af + V2Af. The NFW model is expected 
to be a good fit to the inner parts of the profile, so we use 
the a nalytical expression for th e projected NFW density 
le.g.. IWright fc BrainerdI [200(H) to fit the cluster galaxy 
density profiles within R 200 , establishing the values of 
Us and Tg. We then calculate the logarithmic derivative 
d\og{N)/dlog{R) of the splines for each cluster to test 
for the presence of the density jump feature, comparing 
it to the logarithmic derivative of the NFW fit. 

3.2. Fixed Parameters 

As noted in Section [U we define measures of clus¬ 
ter size such that the mean mass density inside the ra¬ 
dius i?A is p = Apb{z); commonly used values are R 500 
and i? 200 - The other quantity that needs to be spec¬ 
ified is the background density Pb{z)', one choice often 
used in observational work is the critical density Pc{z) = 
3 H{z)^/8ttG, where H{z) is the Hubble constant at red- 
shift z, which is given by = Hq [ 0^(1 + z)^ -I- Ha] 

with Hq = 100 h km/s/Mpc. Another choice is to use 
the mean matter density Pm{z) = pc{z)0,m{z), where 
^m{z) = flm(l + z)^/ [Hm(l + z)^ + Ha] , which is used 
with A = 200. 

As noted above, IDK14I use the mean matter density 
Pm{z) to define the radius R 200 = Rgoom used in the 
fitting formula given by Equation O- However, Fml 
measures R 200 = R 200 C for their sample of clusters using 
the critical density as r eferenc e, so we need to convert 
their measure to that of lDKl4 To do so, we note that a 
given mean density p may be written in two ways: 

p — AcPci^z) = A^p^(z). (8) 

If we specify that p = pnfWj then in the outskirts of the 
cluster (i.e., including near i? 20 o)) we have: 

p(^a)ocT. (9) 
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Fig. 3.— Plots of the projected densities of galaxies with r < 20 in four cluster fields: A1068, A1302, A646, and A1132. Both the 
total galaxy density (black points) and the cluster density (red points) are expected to rise steeply towards the cluster center. The blue 
points show the density after removing cluster members from the total count, with a line drawn to indicate the average value beyond 20 
arcminutes for comparison. The non-cluster galaxy distribution should be roughly flat if we have adequately selected cluster members, as 
for A1068 and A1302. A646 and A1132, on the other hand, show some residual overdensity. 


Combining this relation with Equation ([8)) yields 

A200m _ / 1 y 

For z = 0.1 — 0.3, this implies that 

^200m ~ 1.4 X i?200c. 

Accordingly, for simplicity, we wi ll he nceforth refer to 
^ 200 m as i? 200 j and the values from IRlSi will be converted 
to this measure using Equation (HU. 

Lastly, we need to establish an appropriat e value for 
rim- We note that the original prescription of lDK14l de¬ 
fines Equation © in terms of mass densities p rather 
than number densities n, and their fitting results assume 
a fixed value of Pm{z) = pc{z)flm{z). The translation 
into a number density needs to take into account the im¬ 
pact of the primary selection function by which we obtain 
cluster members, the red sequence method. In the ab¬ 
sence of a cluster, this method would select a population 


( 10 ) 

( 11 ) 


of galaxies that lies in the appropriate region of color- 
magnitude space; in this case, the projected density of 
these galaxies is expected to be roughly constant across 
the field of view. 

Recalling the outer profile term of Equation o, 


rioutir) = rin 


5i?200 


+ 1 


( 12 ) 


we can analytically determine the contribution to the 
surface density of the last, constant term. The surface 
density is the line-of-sight integral. 


N{R) = 2 



r(r)r 


Vr2 - i?2 


dr: 


(13) 


this integral diverges for a constant n{r). However, in 
practice we must truncate this inte gral at some maxi¬ 
mum radius i?max. As noted above. IDK14I use Rmax = 
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10i?vir ~ 9i?200- In that case, 


Nm{R) = ‘2nm [ 

JR 


9-R200 


Vr2 - i?2 


dr, 


Nm{R) = 2nmR200 




(14) 

(15) 


For the scales of interest in our fits, R < 2 i? 2 oo (and even 
a bit beyond), this value is roughly constant, 

Nm « 17.6nmi?2oo- (16) 


The redshift dependence of Nm can be determined by 
applying the red sequence method in test fields that are 
not centered on low redshift clusters. If we select a pop¬ 
ulation of galaxies with this method, then we can con¬ 
struct the projected density in a given radial bin i as 



where A/) is the number of galaxies in the fth bin and 
Da{z) denotes the angular diameter distance. We expect 
that the angular projected density (the term in brackets) 
over the radial range is a roughly constant value, which 
we denote by 77 . Then 


Nm{z) 


V 

Da{zY 


77(1 -I- zY 

D,{zY 


(19) 


where is the comoving distance, which at the small 
redshifts considered here, is given by Dc{z) ~ (c/Hq)z. 
Upon absorbing the factor of c/Hq into 77 , we have: 


iV^(z)=77^i±^. (20) 

z^ 

To obtain the value of 77 , we apply our red sequence cuts 
in random test fields from SDSS (also subject to our ini¬ 
tial magnitude cut of r < 20) and fit iV with the above 
function. This yields 77 « 0.08 /i^ Mpc“^. 

Accordingly, combining Equation (fTCl) with (1^ . we 
find that: 


nm{z, R 200 ) 


4.55 X 10-3 (1 + z)^ 
R 200 z^ 


( 21 ) 


We fix this value individually for each cluster using its 
measured redshift and R 9 .nn (the latter converted as dis¬ 
cussed above) from iRlSl 


4. RESULTS 
4.1. Fitting 

The results of fitting the projected number density pro¬ 
files with Equations m and o are summarized in Ta¬ 
ble [ 1 ] for both the fiducial method and variations, the 
latter of which will be discussed in the next section. 
However, as the DJ model has one additional parame¬ 
ter than the NFW-fOuter model, we then compare the 
two fits using the AlC and BIC. We compute, in each 
case, AlC = 1Cnfw-i-o — ICd j ■ Since lower values of 1C 


are favored, this quantity will be positive if the criteria 
indicate evidence in favor of the D J model. In Table[Tl we 
thus first list in the column, ‘% y^j/ndf < 
the number of galaxy clusters for which the reduced 
is lower for the DJ model and second the num¬ 
ber of galaxy clusters whose fits pass a general qual¬ 
ity cut, including requiring that y^/ndf < 3 for the DJ 
model and a more generous x^/ndf < 5 for the other. 
The next two columns indicate how many clusters have 
AlC = 1Cnfw-i-0 — ICdj > 0 for each criterion, which 
would suggest evidence in favor of the DJ model. The 
last two columns indicates how many clusters have one 
of the AlC > 5, indicating strong evidence in favor of 
this model. 

Our hducial method uses a binning of Ai? = 
0.1 h-^Mpc and employs a magnitude cut of r < 20 . 
We use two selection methods to identify cluster mem¬ 
bers; as discussed in more detail in Section!^ Selection 
A is based on SPSS photometry and Selection B refines 
the first with IR13I spectroscopy. We see that refining 
the cluster member selection yields a higher number of 
clusters for which the DJ model is a better description 
than a simple NFW model with an outer term. How¬ 
ever, the results of individual clusters are consistent be¬ 
tween the selection methods - the seven clusters with the 
strongest evidence in favor of the density jump using Se¬ 
lection A (A655, A1033, A1246, A1437, A1689, A1914, 
and A2034) continue to have AlC > 5 using the second 
selection. These clusters (along with the more marginal 
case A1835, which has only one AlC > 5 using Selection 
A) are shown in Figure|3]for Selection A fits and FigurejS] 
for Selection B; however, using Selection B several other 
clusters also reveal strong (and for some even stronger) 
evidence for a density jump. In both cases, the cluster 
with the strongest evidence for a density jump is A1689. 

Furthermore, the increase in the number of clusters 
with 1C evidence in favor of the DJ model using Selection 
B suggests that it may be benehcial to investigate these 
profiles further with a cleaner cluster member selection, 
with either a sample of galaxies with high-quality pho¬ 
tometric redshift estimates or with a sufficiently dense 
sample of galaxies with secure spectroscopic redshifts. 
Additionally, a few of these clusters (particularly A1246) 
show some residual overdensity at roughly 9 < 5' oi the 
cluster center after the cluster member selection (c.f. Fig¬ 
ure 0 ), which could be indicative either of having missed 
some population of galaxies in our selection or that there 
are some other agglomerations of galaxies along the line 
of sight that contribute to the overdensity. High-quality, 
dense redshift estimates would help resolve this ambigu¬ 
ity. 

We additionally test the fits by making the radial bins 
twice as large (Ai? = 0.2 h-^Mpc). In this case, the 
reduced values are a bit worse for Selection A, as 
indicated in Table [T] We find that the larger binning 
yields the greatest drop in the number of clusters with 
evidence favorable towards the DJ model using the AIC; 
the BIC results, however, are virtually unchanged from 
the earlier case. The discrepancy is likely caused by the 
smaller number of data points over which we fit relative 
to the number of parameters, which makes affects the 
correction term for the AIC. 
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Fig. 4.— Plots of the projected number density profiles for the 8 clusters with the highest AAIC and ABIC values in the fiducial analysis 
using Selection A. Three fitted functions are shown: a base NFW, fitted interior to R200^ NFW+Outer model, given by Equation I® 
and fitted interior to 2R2qo, and the full Density Jump model given by Equation^with an inner NFW profile, fitted interior to 2/?200- 
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Fig. 5.— Plots of the projected number density profiles for the 8 clusters with AIC > 5 using both selection methods, with fits shown 
for the profiles constructed using Selection B. The curves are the same as in Figure [4| 


TABLE 1 

Information Criteria Results 


Method 

% XDj/ndf < Xnfwo/'^'^^ 

AAIC > 0 

ABIC > 0 

AAIC > 5 

ABIC > 5 

Fiducial, Sel. A 

38/56 

18 

22 

8 

7 

Fiducial, Sel. B 

43/56 

32 

32 

16 

20 

AR = 0.2 /i“^Mpc, Sel. A 

34/55 

12 

22 

5 

7 















































































































Fig. 6. — The predicted signature of the density jump is a steep¬ 
ening of slope that can be analyzed using the logarithmic derivative 
of the lDKl^ model for four combinations of Sg and be values. The 
reds hift dependence enters via nm{z, R 200 ) as given by Equation 
123, and we show the result for three values of 2 : in the range of 
the [RQ clusters at fixed F = 4. For the 2 = 0.15 case we also 
show the variation of the profile with F. An NFW profile (dashed) 
is drawn for comparison. All other model parameters are fixed for 
all curves shown. 


4.2. Smoothing Splines 

In addition to fitting the projected galaxy density pro¬ 
files with Equations © and o, we also fit a smooth¬ 
ing cubic spline as described in Section IXTI This latter 
approach has the benefit of being model-independent; 
we use the spline fits to construct a smooth logarithmic 
derivative of the profile to test for the existence of a den¬ 
sity jump. 

First, to provide context for the spline results, in Fig¬ 
ure El we plot the DJ model at fixed T = 4 for various 
redshifts, which determine rim, and for an example fixed 
redshift z = 0.15, we also show the variation with T. The 
bottom panel of each plot shows the corresponding loga¬ 
rithmic derivatives. The effect of the redshift dependence 
is to make the density jump steeper at higher redshift due 
to the lower background rim- At a fixed rim, a smaller 
value of r means a shallower density jump, whose maxi¬ 
mum amplitude slope is attained at larger radii, although 
the values of be and Se, which are free parameters in our 
fits, also contribute to the variation in the amplit ude. 

FigureElcan be compared to Figure 13 of lDKl4 While 
the behavior of the density jump is qualitatively the 
same, we note that in our analysis, the maxi mum sl ope 
amplitudes are smaller than those predicted bv lDKl'j due 


to the elevated background; accordingly, in our case the 
maximum slope amplitudes are governed not only by the 
mass accretion rate T but also the redshift of the cluster. 
However, the density jump location also depends on T. 

Accordingly, we model our spline fits on Figure El We 
show the cubic spline fits for the same subset of clusters 
selected via the information criteria in the fiducial anal¬ 
ysis. Figure [7] shows the spline fits for Selection A, while 
Figure [ 8 ] is its counterpart for Selection B. The spline fits 
do appear to pick up a modest steepening of slope, con¬ 
sistent with the low redshift of these clusters. Since at 
low redshift we can expect to almost exclusively detect 
very large density jumps (T Ri 5), which are likely un¬ 
common, it is perhaps not surprising that the number of 
clusters for which the information criteria suggest strong 
evidence for the jump is fairly low. 

5. DISCUSSION 

We have tested whether cluster galaxy density pro¬ 
files show evidence for a density jump feature near the 
virial radius using two methods: profile fitting and spline 
smoothing. We have examined the evidence in favor of 
the presence of a steep density jump in the galaxy den¬ 
sity profiles of clusters, and also investigated the results 
via spline smoothing. 

There does appear to be some dependence of our re¬ 
sults on cluster redshift. In our fiducial analysis of Sec¬ 
tion [TTI the seven clusters that showed strong evidence 
(both AIC > 5) for the density jump using Selection 
A spanned the redshift range of roughly z = 0.1 — 0.2, 
which omits the higher redshift clusters (although an 
eighth cluster, A1835, with one AIC > 5 value, is at 
z = 0.25). However, employing Selection B enlarges the 
sample of strong evidence clusters, and extends t his to 
about z = 0.1 — 0.25, effectively the entire range ofES 
More massive clusters tend to show stronger evidence for 
a density jump: these same 8 clusters have virial masses 
in the range Mvir = (2 — 11 ) x 10 ^^ while a sig¬ 

nificant fraction of the clusters in the lR13l sample have 
Mvir = (0.3 — 2) X 10^^ h~^MQ (and may also be con¬ 
sidered galaxy groups). This trend continues even when 
considering the additional clusters that pass the crite¬ 
ria using Selection B. This behavior can be expected, as 
higher mass clusters tend to have higher values of T and 
their profiles are likely to be better sampled than lower 
mass systems. 

We see variations based on the cluster member selec¬ 
tion. Our primary selection uses only photometric data, 
but redoing the analysis with the inclusion of the IR13I 
spectroscopy does shift the results; in particular, the re¬ 
fined selection yields a larger sample of clusters that show 
some evidence for the jump. This suggests that addi¬ 
tional data is required to make a firm detection of the 
density jump. This additional data could be in the form 
of dense spectroscopy out to l arge radii of cluster mem¬ 
ber galaxies, building upon the lR13l catalog, which would 
provide the most secure cluster member determination. 
Otherwise, high-quality photometric redshifts could also 
improve upon our estimates. 

Additionally, it is worth noting that in interpreting 
these results it is necessary to keep in mind some of 
the other factors that can impact the jump signature. 
In particular, cluster asphericity and contamination by 
other small groups and clusters along the line of sight 
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Fig. 7. — The top panel shows spline fits for various values of the smoothing parameter S for the clusters in Figured The lower panel 
displays the logarithmic derivatives of the splines. An NFW fit and its logarithmic derivative (dashed line) are shown for comparison. 


can contribute to a diminishing of the jump signal. The 
latter of these can be addressed with redshift estimates. 
For the first, we note that the signature of the jump be¬ 
comes more pronounced and thus easier to detect when 
r is large; however, clusters with large values of F may 
be more likely to have disturbed shapes and substruc¬ 
tures, which when examined in projection could obscure 
the signature of the density jump. 

More generally, it would be useful to compare the re¬ 
sults in Table[T]to simulations. Mock halo catalogs could 
give an indication of the conditions needed for a cluster 
to show a discernible jump and thus the fraction of clus¬ 
ters in which we can expect to find this signature, which 
would provide context for the fractions we find using ob¬ 


servational data. While this is a promising avenue for 
future analyses, such a comparison is beyond the scope 
of this work. 

6. SUMMARY AND CONCLUSIONS 

Using the cluster sample of lR13l and optical data from 
SDSS, we have searched for the signature of a density 
jump feature near t he viria l radius R 200 ) predicted 
by the simulations of lDKl4 Our fiducial analysis selects 
cluster members from SDSS photometric catalogs using 
the red sequence method and photometric redshifts, and 
compares this to a select ion refined by the inclusion of 
the spectroscopy of IR13I . After constructing the radial 
density profiles of the clusters, we fit two models, one 
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Fig. 8.— Same as Figure [7] but for Selection B profiles. 

with a density jump — Equation © — and one without 
— Equation o — and used the Akaike and Bayesian 
Information Criteria (AIC and BIG) to examine the evi¬ 
dence in favor of the density jump model. These criteria 
indicate that, using our fiducial methods, at least 10% 
showed strong evidence (AIC > 5) for the model that 
includes a density jump. The clusters with strong evi¬ 
dence for the density jump tend to have a higher mass 
(-^vir ^ 2 X as expected. The cluster with 

the strongest evidence for the jump in this sample is 
Abell 1689. 

We examined varying some of the parameters in our 
analysis. We find that using a larger bin size yields fewer 
clusters with strong evidence (if we stipulate a strong 


result for both criteria; the AIC appears to be most af¬ 
fected by the binning) for the density jump. However, if 
we require only one of the information criteria to yield 
strong evidence, then we are in excellent agreement with 
the fiducial analysis. 

We additionally tested our results using cubic smooth¬ 
ing splines. The spline analysis appears to indicate some 
steepening of slope in these clus ters, bu t the comparison 
to the theoretical predictions of lDKldl is complicated by 
background and projection effects. In the case of the 
former, a high background can diminish the signal, while 
the latter acknowledges that since clusters can have di¬ 
verse shapes, the density jump is not necessarily localized 
at a single radius, which can lead to modifications of the 
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signal in projection. 

Finally, varying the cluster member selection methods 
does appear to have an effect on the results; in particu¬ 
lar, including the most secure cluster member indicators 
(spectroscopic redshifts) increases the number of clusters 
for which we have strong evidence in favor of a model 
with a density jump. Thus, our conclusions are lim¬ 
ited by the availability of data for the identification of 
galaxy cluster members. While the red sequence method 
provides a means of identifying the majority of cluster 
members, which tend to be red, early type galaxies, it 
does not readily provide a means of rejecting interloping 
red galaxies from higher redshifts or for including the 
smaller population of member galaxies blueward of the 
red sequence. While the spectroscopic observations of 
lR l.1l provide an excellent start, it would be useful to ex¬ 
tend their samples with additional redshift estimates to 
improve the completeness in the cluster outskirts, mit¬ 
igating observational selection effects that can mimic a 
density steepening. 

Accordingly, at present, neither the extant photome¬ 
try nor spectroscopy of these clusters provides adequate 
numbers of confirmed member galaxies out to sufficiently 
large radii to significantly improve upon our methods. 
However, in order to better constrain the density jump 
feature and precisely determine its amplitude and posi¬ 
tion, future efforts will require either secure photometric 
redshifts or additional spectroscopy to refine the cluster 
galaxy density profiles. Future work using mock halo 
catalogs could help pinpoint the conditions under which 
clusters can be expected to show a discernible density 
jump to help guide these observational endeavors. 
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